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Abstract 

We study a renormalized coupling g and mass m in four dimensional 
theory on tori with finite size z = mL. Precise numerical values close to the 
continuum limit are reported for z = 1,2,4, based on Monte Carlo simu- 
lations performed in the equivalent all-order strong coupling reformulation. 
Ordinary renormalized perturbation theory is found to work marginally at 
z = 2 and and to fail at z = 1. By exactly integrating over the constant field 
mode we set up a renormalized expansion in z and compute three nontrivial 
orders. These results reasonably agree with the numerical data at small z. 
In the new expansion, the universal continuum limit exists as expected from 
multiplicative renormalizability. The triviality scenario is corroborated with 
significant precision. 
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1 Introduction 



The standard view |1] , |2] of the primary textbook example of a scalar field theory 
with quartic self-coupling in four space-time dimensions is that it is a trivial the- 
ory. This means that a true continuum limit of the regularized theory inevitably 
leads to a non-interacting Gaussian theory. More precisely, the perturbative renor- 
malization group predicts that, once one is close to this limit, the renormalized 
coupling asymptotically vanishes logarithmically with the cutoff. This weak depen- 
dence keeps 0^ theory useful as a physical theory, because we can simultaneously 
have only small cutoff effects but still have sizeable interaction in an effective 
theory valid over a limited but large range of length scales. In the absence of 
rigorous proofs of triviality in four dimensions, it remains to follow the strategy of 
[2] and to verify in a nonperturbative lattice calculation that the perturbative sce- 
nario actually applies close to the continuum limit, and thus to exclude the logical 
possibility that it is simply irrelevant. We re-address this old subject because nu- 
merical progress [3] has made it possible to render such tests more stringent with 
only moderate resources. Although the situation is complicated by the coupling to 
other fields, triviality is also the standard scenario for the Higgs field in the Stan- 
dard Model. The subtle nature of logarithmic effects makes a numerical check of 
the triviality picture rather non-trivial in spite of the simple field structure of 0"^ 
theory. A number of Monte Carlo simulations have been conducted of which [1] 
is one example where the then new cluster algorithm was used in the symmetric 
phase. These as well as our new simulations took place in the Ising limit of (p^ 
theory, which amounts to infinitely strong bare coupling, where we expect to find 
the strongest possible interaction. 

In [;3j a novel approach to simulate this model was presented which is based 
on the simulation of Aizenman's random current representation [5]. Its numerical 
efficiency is closely related to the fact that this equivalent reformulation enabled 
Aizenman to derive bounds that prove triviality in D > 4 dimensions. For both 
applications it is essential that the connected four point correlation that enters 
into the renormalized coupling, can be calculated without the need of performing 
cancellations and with the Lebowitz inequality |i6jj being manifest. The elimination 
of the otherwise unavoidable loss of significance in the Monte Carlo allows to rather 
easily compute the renormalized coupling with per mille precision close to the 
continuum limit. An additional decisive bonus is the practically complete absence 
of critical slowing down. 

The triviality conjecture in D = 4 is investigated here in a finite volume renor- 
malization scheme where the continuum limit is approached on a volume that 
remains finite in units of the physical renormalized length scale. This has the 
same advantage here as in asymptotically free theories [7], [8], namely that the 
mastering of a large ratio of physical scales is avoided and one can thus probe the 
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universal continuum limit much more closely . 

In this paper we employ the numerical technique described in [3] without 
changes and therefore here leave out all details on this. The data with z = 2, 
already published in |3], are extended by including a lattice of size 64^ and by 
simulating two more series of lattices at z = 1 and z = 4. In addition we report 
on a computation in lattice perturbation theory that we employ to analyze the 
data. This calculation goes beyond [2j by allowing for finite z values and by com- 
puting the lattice artifacts of the Callan Symanzik /3-function up to two loops. 
By combining with [2| we also obtain the three loop term for various finite z. By 
comparing with our data we will conclude that perturbation theory works well for 
large z and fails for small z. 

A plausible reason for the failure of perturbation theory on a small torus is 
that the constant zero momentum mode may receive too little Gaussian damping 
to justify its perturbative treatment [9]. We therefore embark on an alternative 
approximation scheme where we treat this one mode exactly while maintaining 
the perturbative expansion for all others. It turns out that this type of expansion 
rearranges itself under renormalization to an expansion in powers of for arbitrary 
values of g/z'^. The new expansion is found to be applicable and accurate at small 
z. 

We set up our renormalization scheme in sect. 2. Our perturbative calculation 
is described in sect. 3 and app. A and numerical results are discussed in sect. 4. 
The new small volume expansion is outlined in sect. 5 and app. B and we conclude 
in sect. 6. 

2 Finite size renormalization scheme 

We define 0^ theory on a four dimensional periodic lattice of extent L in all direc- 
tions by the standard Euclidean action 



Here = [(j){x + afi) — (^[x)\la is the standard nearest neighbor forward 

derivative and /i is a unit vector in the /i direction. From here on we mostly use 
lattice units and leave out powers of a = 1 except for some formulas where we find 
that clarity is gained by restoring explicit factors of a. A completely equivalent 
form of the above action on the lattice is given by 



S = a 





S = -2kJ2 s{x)s{x + fi) + + A(s(x)2 - 1)2]. 



(2) 
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The relation between the two parameterizations is given by 

^o = 6A/k^ = (1 -2A)//t-8, 0=^2^5. (3) 

In the latter form it becomes manifest that for A — )■ oo the integrations over the 
spin field s{x) reduce to Ising sums over s(x) = ±1. 

It has become standard [2| to take the continuum limit along vertical lines in 
the {k, A) plane by sending n to its critical value Kc(A) at fixed A. If for the infinite 
volume this limit is taken from below the theory is in the symmetric phase. We 
adhere to this 'coordinate choice' although alternative procedures are conceivable, 
of course, without however changing the set of continuum theories that can be 
reached. 

It is now a completely well-defined procedure to perturbatively compute, start- 
ing from ([T]) and within the regularized theory, the (bare) effective actior0 as a 
sum over connected one particle irreducible graphs with propagators for external 
lines canceled. At first we are interested in the 2-point vertex function T^'^\p, —p). 
It is related to the susceptibility measured in simulations by 

X2 = 5^(^(0).(x))= [-2«:r(2)(0,0)]-\ (4) 

X 

For our definition of a renormalized mass m a non-zero momentum is required for 
which we take the minimal one 

= (27r/L, 0,0,0). (5) 
Then the 'second moment' definition of m follows from the universal ratio 



(6) 



with = 2sin(p^/2). In the following we define the finite size continuum limit 
by holding constant the combination 

z = mL. (7) 

Thus for some choice of z (and A) we consider sequences of lattices with growing 
L = L/a — i- oo where for each L we adjust k such that z has the desired value. We 
may also say that in this way we have defined a family of renormalization schemes, 
one for each value of z. 



""^The sign of F is not uniform in the literature. We follow the convention which leads to 
r(2)(0,0) = -m2 + 0(5o). 
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The renormalized coupling is now an output observable in this procedure which 
is given by another universal ratio 
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with the connected four point susceptibihty 



X, = J2{siO)six)siy)s{z))-3L\l 



1 rw (0,0, 0,0) 
(2/€)2 [r(2)(o,o)]4 



(9) 



or, equivalently, the 4-point vertex function T^^\ It vanishes for a Gaussian the- 
ory and is hence a key quantity in connection with the triviality conjecture. An 
important property of our coupling is that it is rigorously bounded in the range 



This bound, based on [5j, is manifestly visible in eq. (20) of [Sj where an observable 
with values in {0, 1} is averaged with a positive weight to yield g/{2z^). 

The renormalization scheme in [2] is defined in an infinite volume and uses the 
'zero momentum' definition for the renormalized mass m = mji. We make contact 
with this scheme by taking within our family of schemes the simultaneous limit 
2; — cxD, L — )■ oo at fixed rriR = z/L. In this limit pl/m^ ~ {211/ zY becomes 
arbitrarily small and ([6]) goes over into (2.12) of [2]. Note that this is not a 
continuum limit as long as = airtji is finite, in contrast toL = L/a— T-ooat 
fixed z. 

We are now in a position to define the Callan Symanzik /3-function in our 
scheme by 



As L is integer, the derivative must be approximated by a sufficiently accurate 
finite difference formula, see sect. 13.11 for further details. 

There is a family of 'curves' in the ((?, a/ L) plane, one for each value A G [0, 00), 
on which and other observables become defined. We write 'curve' because 
with only integer L/a it is actually a sequence of discrete points. If the triviality 
conjecture holds, then all these curves end in the point (0,0). In figure [1] we see 
a sketch of the expected shape of the domains for z = 4 and z = 2 based on our 
nonperturbative data together with continuum extrapolations to be discussed in 
section |H The respective domains are the areas under the curves if we assume 
that smaller A yield smaller g. An analogous curve for z = 1 would be roughly 
around g ^ 1.5 before diving down to zero. 



^ g/z^ ^ 2. 



(10) 



Pz{g,a/L) = -L—g\x,z- 



(11) 
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Figure 1: Upper boundaries of the domains of renormalized coupling versus cutoff 
at z = 2,4. Symbols are Monte Carlo results and 1-loop perturbation theory 
furnishes the dotted extrapolations. 



At tree level of perturbation theory we trivially find g = go and m = mg. In 
terms of the natural variables in f llip this reads however 

^ = ^^^(8 + .V^Y (12) 

which leads to 

4^2/7-2 

p. = ',,T^ 9 + 0(^7^) (13) 
8 + /L 

as already noticed in [2]. To avoid this tree level lattice artefact in the /9-function 
we change the definition of the renormalized coupling [3] by a term that vanishes 
quadratically with the cutoff into 

~g = gr\z/L) (14) 



with 



The /9-function for this renormalized coupling 

'dL 



P,{~g,a/L) = -L^g\x,. (16) 
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vanishes at tree level. 



3 Lattice perturbation theory up to two loop or- 
der 

3.1 Artifacts of the one and two loop beta function 

In this section we discuss the asymptotic expansion of (3^ in powers of g 

/3,aa/L) = 5^6«(a/L)^'+^ (17) 

i>i 

or the completely analogous formula without the tildes [which must include / = 
however]. The (perturbative) renormalizability implies the finiteness of all limits 
\ima/L^obz\a/ L). The one and two loop terms are scheme independent, which 
here means independent of z, and have the universal values 



hm W\a/L) = = At, lim W\a/L) = ^ = (1^ 



In [2] the three loop result 



is given for the infinite volume scheme. 

In appendix |A] we derive the coefficients of the following expansion 

~g=~gQ+ Piiz, L)gl + p2(^, L)gl + O(^o') (20) 

with 

9o = Qor^imo)- (21) 

We have performed the necessary Feynman diagram sums up to L = 100. The only 
technicality that is perhaps worth mentioning here is that by a judicious use of 
both momentum and position space propagators and the fast Fourier transform we 
compute all two loop diagrams by performing no more than 0{L'^\nL) operations 
at each L. More details are given in app. |Al 

Because of the relation g^ = 384A/(1 — 2A)^ we may keep constant in (ITB]) 
instead of A. Then we obtain 

WHl-') = -L-^Mz,L), hf\L-') = -L-^[p,{z,L)-Mz,Ln (22) 
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For several z we have checked that in the Symanzik expansion in terms of In' LL~^^ 
the exact values fe^^'^) of (fT8|) emerge as leading terms (/ = 1, n = 0) with significant 
precision. With this verified we cancel them and form the deviations 



(23) 



(24) 



d_ 

bW^dL 

m dL 

We use the four point formula for the derivative which contributes errors of order 
to di^^ which itself is expected of size 0(L~^). Results are shown in Fig. [2j 
The dots are the corresponding lattice sums while the lines are fitted third degree 
polynomials in to interpolate and represent the data. 
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Figure 2: Deviation from the universal values of the 1 and 2- loop /3- function at 
finite lattice spacing. 



3.2 Three loop beta function 

We now want to relate the renormalized couplings g and g' referring to schemes 
defined with two values z and z' that we freeze for a while. To appreciate the 
relation we first give an 'operational' description. We first choose a value L = 
L/a ^ 1. Then there are (assumed to be) bare parameters A, k, in the scaling region 
that lead to renormalized values z = mL and g. Now with the bare parameters 
unchanged we change L — )■ L' ^ 1 until z' is found. In this way g' and L'/L 
become functions of g up to scaling violations of order that we neglect in this 
subsection. With these relations known the /3 function transforms as 
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In the end we shall be interested in the limit 2;' — ?• oo with g' — >■ to make contact 
with the infinite volume scheme for which we know the 3 loop /3 function (|T9|) . 

The relations just described are now analyzed in perturbation theory. Then 
there is an expansion 

g = g' + P^iz, z')g'' + P^{z, z')g'^ + 0{g''). (26) 

From this form it follows that the first two expansion coefficients of the /3-functions 
are independent of z while at three loops the connection is 

If - bf = 6(1) (P2 - Pi') - Pib'-'^ . (27) 

Using (155]) to leading order we first derive 

L'/L = z'/z[l + Q,{z,z')g' + 0{g'')] (28) 

with 

Q,{z,z') = l^[q,[z\Lz'/z)-q,{z,L)] + 0{L-^). (29) 

Then we find 

P,{z,z')=p,-p\ + 0{L-^), (30) 
P2 = P2 - P2 - - Pi) + h^^'^Qi + 0(L-2). (31) 

Here Pi^p'i denote Pi{z,L) and Pi{z' , Lz' / z) respectively. Arguments L' have been 
eliminated with fl28p and use was made of Pi{z', L') = —b^^^ InL' + 0{L'^). 

We have defined now a number of expansion coefficients that emerge as finite 
continuum limits as L — )■ 00. The Symanzik expansions of the corresponding 
combinations of lattice Feynman diagrams have been analyzed by the method 
given in appendix D of [lOj. Results are collected in Tab. [H We see that the 
3-loop coefficient bf rises very steeply for -2^3. The negative entries in the 
last column imply that as one lowers z starting from z = 16 (which is effectively 
infinite) to smaller values, the behavior of bf is not completely monotonic. One 
might be tempted to think of numerical inaccuracies here, but as far as we can tell, 
the negative sign seems to be significant. The table allows to change z in cycles 
for which we find consistency. A rough picture is that bl is constant above z = 4 
and rises by a factor 15 for 2; = 4 — )■ 2 and by another factor 270 for z = 2 — > 1. 
The typical couplings will at the same time be seen below to diminish by factors of 
roughly 1/3 and 1/7. This means that the 3-loop contribution overall rises steeply 
compared to the 2-loop term with its ^-independent coefficient. In this way the 
perturbative series indicates its breakdown for small z. 



9 



z' 


Qx X 10^ 


Pi X 10^ 


Pi X 10^ 


{hT - h^/) X (47r)*5 


1, 2 


210.06698 


-1394.9977 


3437.304 


110315.00 


2, 4 


9.7211061 


-80.361921 


12.52668 


382.0849 


4, 8 


0.2310329 


-2.952908 


0.04461 


0.0425 


8, 16 


0.0010522 


-0.024524 


0.00024 


-0.0037 


2, 3 


8.6999022 


-70.199600 


10.75077 


372.7864 


3, 6 


1.2399938 


-12.901805 


0.44644 


9.4008 


6, 12 


0.0132854 


-0.237614 


0.00203 


-0.0651 


4, 6 


0.2187899 


-2.739483 


0.04164 


0.1022 


8, 12 


0.0010424 


-0.024189 


0.00022 


-0.0053 



Table 1: Expansion coefficients relating renormalization schemes defined by z and 
z' respectively. Numerical errors are beyond the digits quoted. 



4 Analysis of precise numerical data 

As mentioned in the introduction, our simulations here follow in all details those 
described in |3j. In particular, for each set of parameters we have generated a 
statistics of 10® iterations resulting in per mille errors. Due to the efficiency of the 
method this could be done in a short time on a few up-to-date PCs. 

In the next subsections we list data and plot couplings versus cutoff for z = 
4, 2, 1. The tables contain the quantity X that is related to the coupling by 

g = 2z^X. (32) 

The evolutions of the coupling are compared with curves obtained by integrating 
the renormalization group equation 

L^ = -~M~9,0) (33) 

starting from the 'measured' coupling at the largest L/a = 64. The (solid) curves 
labelled with a loop order refer to the perturbative expansion of (3^. We have 
also looked at the curves including the known one and two loop cutoff effects, 
i.e. with PzigA/L) on the right hand side. We find however that these small 
corrections do not systematically improve the picture. Sometimes they go in the 
right and sometimes in the wrong direction. We below offer a possible explanation 
for this and do not include these curves in the plots. For the case z = 2 they 
are visualized however (without L = 64) in Fig. 2 of |3]. Instead we here include 
additional curves (dashed) labelled LO, NLO, NNLO. They refer to the leading, 
next-to-leading and next-to-next-to-leading orders of the small volume expansion 
that is explained in detail in section EJ 
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4.1 z = 4 



L/a 


2k 


z 


X 


dX/dz 


X{z = 4) 


8 


0.141976 


4.0025(22) 


0.07684(29) 


-0.0556(5) 


0.07698(27) 


10 


0.144491 


3.9956(21) 


0.06783(27) 


-0.0518(5) 


0.06760(25) 


12 


0.145933 


3.9991(20) 


0.06173(26) 


-0.0489(5) 


0.06168(24) 


16 


0.147481 


4.0045(19) 


0.05419(24) 


-0.0441(5) 


0.05439(23) 


22 


0.148454 


4.0004(18) 


0.04690(22) 


-0.0387(5) 


0.04691(21) 


32 


0.1490781 


4.0002(17) 


0.04120(21) 


-0.0360(5) 


0.04121(20) 


64 


0.1495244 


4.0009(15) 


0.03314(18) 


-0.0292(7) 


0.03316(18) 



Table 2: Monte Carlo results for z = 4. 

Data for z = 4 are tabulated in Tab. |2] and plotted in Fig. |3l At this z value 
perturbation theory works as one expects: the successive terms alternate around 
the nonperturbative answer and get closer to it. The change in the coupling as 
the cutoff changes from L/a = 64 to L/a = 32 is reproduced with 1.1% error in 
the 3-loop approximation (2.4% and 3.9% for 2- and 1-loop). In the small volume 
expansion LO and NLO fall almost on top of each other and are close to 1-loop 
while NNLO is not far from the 2-loop result. 



4.2 z = 2 



L/a 


2k 


z 


X 


dX/dz 


X{z = 2) 


8 


0.148320 


1.9981(27) 


0.39235(96) 


-0.3200(14) 


0.39175(63) 


10 


0.148748 


1.9949(26) 


0.37256(92) 


-0.3193(14) 


0.37093(62) 


12 


0.148996 


1.9992(26) 


0.35493(91) 


-0.3165(15) 


0.35469(60) 


16 


0.149270 


1.9988(25) 


0.33161(91) 


-0.3129(16) 


0.33125(58) 


22 


0.149449 


2.0085(24) 


0.30831(86) 


-0.3030(16) 


0.31088(57) 


32 


0.149571 


1.9956(24) 


0.29028(83) 


-0.2993(20) 


0.28896(55) 


64 


0.1496564 


1.9893(22) 


0.26016(78) 


-0.2887(27) 


0.25706(51) 



Table 3: Monte Carlo results for z = 2. 

Data for z = 2 are tabulated in Tab. Eland plotted in Fig. |H For z = 2 the one 
loop result happens to be exact within errors (0.2%), while 2- and 3-loop are off by 
1.4% and 3.9% respectively for the same benchmark as discussed in the previous 
subsection. As already concluded in [3] this is a limiting case for perturbation 
theory as an asymptotic expansion, where higher loops do not help at all any 
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Figure 3: Evolution of the the couphng g with the cutoff L/a for z = A. 



more. The precision of the small volume expansion is 2.6%, 2.4%, 1.0% for LO, 
NLO, NNLO respectively. 



4.3 z = l 

Data for 2; = 1 are tabulated in Tab. |4]and plotted in Fig. |5l At 2; = 1 perturbation 
theory is not useful anymore. The coupling itself (but not g/ z'^ on its natural scale) 
has small values and therefore the 2-loop contribution, whose coefficient cannot 
depend on z, is tiny. The very large 3-loop term knocks the approximation far 
to the other side. The small volume expansion exhibits deviations of 0.7%, 0.5%, 
0.3% for LO, NLO, NNLO. 

4.4 Lattice artifacts and perturbation theory 

In this subsection we offer a short discussion whether or not we expect the lattice 
artifacts to be well represented by perturbation theory. This may be questionable 
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l/ln(L/a) 

Figure 4: Evolution of the the couphng g with the cutoff L/a for z = 2. 



even in the regime where the continuum behavior is reproduced. 

We want to analyze closer the perturbative procedure that we have applied, 
which is however the conventional one. To derive the expansion coefficients f|T7|l 
we consider qq as an arbitrarily small expansion parameter and correspondingly 
expand around = in ([T]). Then we re-expand the truncated series in the 
renormalized coupling g ot g. This series is then held against lattice data produced 
with A or go such that the bare lattice action has a pair of nonzero constant 
minima. This is so in the most extreme way for the Ising limit. Renormalized 
perturbation theory is expected to describe the universal physics in spite of the 
apparent contradiction as long as successive terms in the g expansion (at small 
a/L) look reasonably 'convergent'. A possible explanation is as follows. We may 
simultaneously consider a lattice that is coarser but still in the scaling region 
tuned to the same (matched) renormalized parameters. On this lattice the bare 
go will be smaller and the expansion may also be naively justified. Up to differing 
small scaling violations the universal physics will then agree on both lattices. 
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L/a 


2fi; 


z 


X 


dX/dz 


X{z = 1) 


8 


0.151670 


1.0020(25) 


0.79953(84) 


-0.4444(30) 


0.80040(85) 


10 


0.150900 


0.9997(24) 


0.78776(88) 


-0.4641(30) 


0.78764(83) 


12 


0.150498 


0.9999(24) 


0.77930(92) 


-0.475(30) 


0.77926(82) 


16 


0.150118 


1.0032(22) 


0.76121(96) 


-0.5074(31) 


0.76281(80) 


22 


0.149907 


0.9961(22) 


0.75058(101) 


-0.5213(32) 


0.74854(78) 


32 


0.149787 


1.0032(21) 


0.72936(106) 


-0.5436(50) 


0.73110(75) 


64 


0.1497143 


0.9993(20) 


0.70179(110) 


-0.5797(48) 


0.70141(71) 



Table 4: Monte Carlo results for z = 1. 

One could consider the extension of this matching to the leading cutoff effects 
in the spirit of the Symanzik effective action. But then one would need further 
improvement terms on the coarser lattice to reproduce the cutoff effects of the 
finer one. As an alternative qualitative argument we could imagine a block-spin 
coarsening of the original lattice. Then we would expect to arrive at an effective 
action on the coarse lattice consisting of the relevant terms of ([T]) with small Qq and 
additional contributions corresponding to Symanzik terms. Therefore the relation 
between cutoff effects computed in perturbation theory with the standard action 
(II]) as sketched above and those of the nonperturbative results may not be so 
straightforward in our opinion. 



5 Small volume perturbation theory 

In this section we introduce a modified perturbative expansion which will turn out 
to lead to a systematic expansion in the finite size scaling variable z defined in ([7]) . 

In our perturbative calculation the Gaussian damping of the low momentum 
modes on the torus is controlled by u"^ = + {2TT/L)'^fP with n having small 
integer components. For small z with itlq = z/L + 0{go) the n = mode receives 
little damping and ordinary perturbation theory may be a bad starting point [H]. 
We therefore split up the lattice field 

^{x) = ^jj^^ + r]{x), y2v{x) = 0. (34) 

Lgo 

We decompose the action ([1]) into three parts, S = Sq^ + S'o,.^ + Si, 

So^i = + (35) 
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l/ln(L/a) 

Figure 5: Evolution of the the couphng g with the cutoff L/a for z = 1. 



and 



with 



3/4 

+ } (37) 



° (3^ 



/go 

For every cutoff L/a we consider now ZqjQq as bare input and want to compute 
from them z and g. The idea is to treat 5*0,. exactly and 5*1 as a perturbation, 
i.e. qq is small at finite Zq. In the end the propagator implied by (136|) will also be 

1/2 

expanded in Qq . Alternatively it could be included as a two 1] vertex in 5*1 from 
the beginning which we found less efficient however. 

The second moment renormalized mass definition ([6]) implies 



z 



2 



I + Z^Kl (02) 



^0 ,79\ (39) 



15 



with 



and 



3L2 



(40) 



(41) 



By doing perturbation theory in we first produce the coefficients of the ex- 
pansions 



my 



^CnizlL)gQ 



n/2 



n>0 



and from (139 



n/2 

■ 



The P function is computed as 



For the derivative we use 



d 
'dL 



-L 



d 



dg_ 
dL 



1 2,30- 



d 



12,30 -Qj^ '^^^"'^^5^2 



(42) 



(43) 



(44) 



(45) 



where all partial derivatives on the right hand side are now taken with respect to 
the set {zo,L,gQ) and 



dz^ 
'dL 



dz^ 
dzl 



(46) 



has been introduced. 

If we invert (H3|) to express ^/go as a series in z"^ we obtain 



g/z^ = ^en{zo,L)i 



,2n 



from fl42D and we arrive at the intermediate form 



/9. = ^'$^/n(4^)^ 



,2n 



(47) 



(48) 



Below we now refer to results whose derivation is described in some detail in 
appendix [Bl 

The ffist coefficients Cq, ei are given by 



Co — -J 5- — -J 5 



(49) 
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and 

ei = {^l''zl + (3 -.4 _ 2) _ i8^^2 + 12} = ^^Ci^^e^ 

with the second moment computed with (153]) 

and the constant Ci defined in fllOOp . For given g, z we solve for Zq in 

g/z^ = eo(^o) 



(50) 



(51) 



(52) 



which gives Zq the status of a renormalized parameter. One finds for example 
Zq ~ 0.303 (/i ~ 1.25) for g = 10,z = 2. The bare quantity Zq is then written as a 
power series in by solving fll7|) 



^0 + X1^"'^^0'^)^ 



.2n 



The lowest order terms are 



(53) 



h3 



-ei 



9zh 



62 + h] 



dei 1 (9^eo 



dzl 



+ - 



2 9(zo^)2 





'deo' 




[dzl\ 



^^0 



2 9(^2)2 



6 9(^2)3 



d{zl 



dzl 



(54) 



(55) 



^^°^56) 



dzl\ 



with all quantities on the right hand sides taken at Zq. To finally obtain the /3 
function in the form which is expected to possess a continuum limit we eliminate 



from fj48l) to obtain 



(57) 



where the -B^ are functions of g/z'^ (via Zq) only after neglecting lattice artifacts 
proportional to L^^. This also is the reason why there is no 0{z^) contribution, 
since 

/o = 0(L-2) (58) 

1/2 

holds. This can be understood by noting that the leading order g^ correction is 
essentially a mass renormalization, the term proportional to (j? in fl99p . This leads 
to relations 

1 „ dcQ 



Cl 



di 



1 ddp /i^ + 2zliJi - 2 

2 2ifL/i2 



(59) 
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and also to the last equality in ( 150|) . With these identities the cancellation that 
leaves only cutoff effects for /o can be shown. The first three Bk are 

5i = fi + hM, (60) 



Bo 



dzl 

(62) 

where the limits L — )• oo are understood on the right hand sides. It turns out that 
these combinations are such, that all divergences cancel and only finite universal 
results enter into the continuum limits of the B^, see app. [B]for more details. 
Our first result is 

^ 9 ^%^ + 2/i^(3^o^-2)-18/i^g + 12 

For large Zq the fluctuations of become Gaussian and fi may be expanded as in 
f lllSp . We then find the leading large behavior 

eo~4' ^i-tI^4' ^'^i-M' i^o^oo) (64) 
Zq 167r"= 

such that in this limit we recover the perturbative result with the one loop coeffi- 
cient flTSl) . More generally we may run Zq through some range and parametrically 
generate the graph of Bi versus g/z^ as shown in Fig. |6l 
In the next order we find 

B2 = -^Bl (65) 

At the moment we see this form on the basis of our series expansion for B2 and 
cannot give a simple reason for this simple relation. In the perturbative limit 
we thus obtain B2Z^^ ~ —{4:/9)b1g^/z^ which is a contribution that in ordinary 
perturbation theory could only come at the 3-loop level. 
Finally we find 



167r2 



-i?3 + — /i^ 



47r2 9 



(66) 



with 



^3 = (^2^6/i^^-6)2 + + 4/.6(204 - 1) + 4^5^2(24^4 _ ^3) 

-16/(9z^ + 34 + 4) + A8fi^z^{154 - 7) - 228;u2(54 - 1) 

+1296/iz^ -432}. (67) 
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Also these factorizations have been found by inspection only. The perturbative 
limit 

i?3--4 (^o^oo) (68) 

^0 

yields B^z^"^ ~ 1)29^ with the term with R3 not contributing to the leading order. 
We thus recover the 2-loop term f|T8l) . 

The running with the leading (LO, Bi), next-to-leading (NLO, -81,2) and next- 
to-next-to-leading (NNLO, -81,2,3) order /3 function in the small volume expansion 
is shown by the dashed lines in Figs. |3] , H] and |5l 



6 Conclusions 

We have reported and discussed data for a reasonably defined renormalized cou- 
pling g in (j)^ theory in finite size systems with z = 1, 2, 4 on lattices up to 64^. All 
simulations took place at infinite bare coupling. The dependence of g on the cutoff 
a/L at z = A was found to be well described by perturbation theory with the first 
three loop orders becoming successively more accurate. Thus there is all reason to 
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trust the perturbative continuation of these curves all the way to vanishing lattice 
spacing a which implies triviality. At 2; = 2 on the other hand it was already found 
in [3] that a (presumably accidental) excellent agreement of the nonperturbative 
evolution with the one loop approximation holds, but that it only deteriorates 
when further terms are included. We now have added the finding that aX z = 1 
the standard perturbative description fails completely. 

The suspicion that the constant mode on the torus is the culprit responsible 
for the failure has led us to work out an expansion where this one mode is treated 
nonperturbatively. Upon renormalization this has led us to an expansion of the /3 
function in powers of z itself for finite values oi g/ z^ for which we rigorously know 
that their domain is contained in the interval [0,2]. We construct completely 
explicitly the first three orders of this expansion, neglecting only order 2;^^ and 
higher, and find it to work very well at our smaller z values. Since this approxima- 
tion also yields a positive (3 function, the coupling evolution may be continued to 
vanishing a and triviality follows once more. The intermediate value z = 2 seems 
somewhat problematic for both expansions and one loop truncated perturbation 
theory happens to be better than the z-expansion in this difficult range which was 
unerringly chosen as the first application of the new simulation technique in [3]. 

We end on a more general remark. Many structural properties of quantum 
fields theory like renormalizability and the existence of the continuum limit are 
(in favorable cases) based on proofs to all orders of perturbation theory. It is 
not so often emphasized, that for nonperturbative calculations on the lattice one 
has to assume that this holds true also beyond perturbation theory, although we 
know that the latter can become numerically completely irrelevant, even if one has 
command over high orders. In our small volume expansion which may perhaps be 
seen as capturing a glimpse of nonperturbative behavior in an analytic treatment, 
we could confirm these standard assumptions. 

Acknowledgments. We thank Jean Zinn- Justin for discussions in the early stage 
of this project. U. W. has enjoyed hospitality and support by the Max Planck 
(Werner Heisenberg) Institut in Miinchen. Financial support of the DFG via SFB 
transregio 9 is acknowledged. 

A Perturbative expansion 

We here report details on the computation of pi 2{z,L) entering in section 121 In 
terms of IPI vertex functions our renormalized parameters read 



Z 



-1 



[r(2)(o,o)-r(2)(p,,-p,)]/p2, 
^2 ^ _^r(2)(o,o), 



(69) 



(70) 
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^ = -z2r(^)(o,o,o,o) (71) 

with the wave function renormahzation factor Z and all quantities at finite L. 
Standard bare perturbation theory gives 

r(2)(p, -p) = -ml -f- gJ-Ji + 9l]jiH^ + gl\j2{mlp) + 0{gl), (72) 



r(^)(0, 0, 0, 0) = -go + gf-H, - gl (^H^, + ^-H^, + ^-H,}j + 0{gl). 



(73) 



The capital letters stand for the usual Feynman diagrams for the two and four 
point functions up to two loops. With all external lines at zero momentum except 
in J2(mQ,p) they are regarded as functions of rriQ (and L) at this stage. Their 
actual evaluation proceeds via the following sequence of steps, 

GW = ;^E<'*"G(p), (75) 

P 

Ji(mg) = G(0), (76) 
^i(^o) - YP{P)? = E[^(^)]' (Plancherel). (77) 

p X 

The momentum sums run over all values = 27rn^/L, = 0, 1, 2, . . . , L — 1. 
Next we compute 

G-{p) = J2[G{x)re-'P\ n = 2,3, (78) 



and use it in 



H2a{ml) = ^^^G\p)[G{p)r, (79) 



p 

H,,{ml) = [G2(0)]2 = H,{mlf, (80) 

^2c(m^) = Jl-lE[G(p)]^ (81) 
p 

Umlp) = G\p). (82) 

The Fourier transformations are performed as FFT on one coordinate direction 
after another, schematically like 

G{p) = F{po,PuP2:P3) F'{xq,Pi,P2,Pz) F" {xq, Xi,P2,Pz) . . . . ^ G{x) 

(83) 
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at a total computational complexity of DL'^ In L only. 

With these expressions we can write (omitting the remainders. . . . + 0((7g)) 

^ = 1 + - (84) 



and 



2 2 

.2 _ ™2 ™2 _ 9o J , 9o J Tj , 9o 



3„ 3„ 3„ 1 



Am^ = - = -| + '-fj,H, + | [(1 + m^/p^) - {ml/pDMmlp^)] 

(85) 

^7 = ^70 - ^o^^i + ^0 (^3i/2a + ^^26 + ^^2c + ^[-^2(^2,^0,) - Ja]^ • (86) 

In these expressions the mass is still mo in all diagrams. In order to obtain g as 
a function of go and we have to combine the last two lines to eliminate ml on 
the right hand sides. To the order considered and using 

dJi dHi 



-^1' "^13 — 2 ^ ~2if2c (87) 



we arrive at 
with 



Am^ = qi{z, L)go + q2{z, L)gl 
1 



and then 
with 



gi(^,L) = --Ji(m2) (89) 

g2(^,i^) = ^[(l + mVp^)J2(m^0)-(mVpDJ2(m^p,)], (90) 

^ = ^o + Pi(2;,L)^o +P2(2;,L)^i^ (91) 

Pi(^,L) = -^ifi(m2), (92) 

P2(^,i:) = ^H2a{m^) + -H2b{m^) + ^[J2{m\p,) ~ J2{m\Q)]. (93) 

4 3p^ 

We now finally change from the expansion of g in powers of g^ to the one of g 
in go of fl2U]) with the corresponding coefficients. Some straightforward steps lead 
to 

gi = r"^gi, q2 = r~^q2 + —ll (94) 

and 

Pi = r"V + ^^i, P2 = r"^p2 + — Pigi + ^ (^'?2 + — gi j (95) 
with the mass m in all arguments here. 
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B Small z expansion 

By straightforward manipulations one can show that 

(^(0))o,0 

holds with a polynomial in ^ deriving from 

P(0) = Po(0)(e-^Oo,r, (97) 

1/2 

truncated at the desired order in g^^ . The subscripts of the averages refer to the 
parts of the action f l5^ . f l5B]) used. To construct P we perform Wick contractions 
with the ri propagator 



Ao(p) 



p2 + ^;/2(2^2 + ^2)/(2L2) 



^o(^) = ilE^^''^o(P)' (98) 



1/2 

which will be expanded in in the end. The factor Pq derives from the Gaussian 
integral over e""^'^''' and is given by 



(99) 



with 



(100) 



The dependence of Ck on L will be discussed in more detail below. By computing 
connected graphs we obtain in addition 



ln[P/Po] = -^X2 + ^02L2A3(O) + | L^A^(0) + 3X2Ag(0) ^Q{gT) (101) 



with the same notation as in the (!78|) and the short hand 

fe-i 



X = L2Ao(0) = 5^Cfc 



1/2 

-^(2^0^ + 0^) 



(102) 



To compute the mass in a similar fashion we expand 

(p^A(p,)e-^Oo,, = P(0) X g(0) 



(103) 
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with Q built from connected diagrams with two external lines p — ±p*, 

3/2 _ ^2 



Q 



with 



n>0 



1/2 

^° (2.0^ + 0^) 



2K. 



evaluated up to the required order. We introduce additional constants for 



l/2/r,^2 , 12^■nl 



and 



L''-'A^M = Dl + g]!\2zl + 0^)^*' + 0(^o) 
They are given by 

k ^ 



and 



1 1 



(104) 
(105) 

(106) 
(107) 

(108) 

(109) 



and corresponding formulas for and D"^ with replacing zero for the total 
momentum. Wc now discuss the behavior of these constants as far as they enter 
into our computation. 

We first note that there are a number of universal logarithmic divergences that 
enter into our final expansion coefficients Bk, like 



^ +0{L-')^L^DI. 



(110) 



$7r^ ■ ■ dL 
The constants and are quadratically divergent, but their difference obeys 



Another combination that was found to occur in our final result is 

1 



d_ 
dL 



3 



647r4 



+ 0{L-'). 



(112) 
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All other constants are non-universal (reflect the hypercubic lattice and the dis- 
cretization) and are either divergent {Ci, D4, D3 + D^) or finite. They all drop 
out in our finite universal results for the (3 function. The behavior claimed above 
was checked by computing the constants numerically in the same way as reported 
in the previous appendix. However, the universal coefficients, which are the only 
feature needed here, should also be computable more easily in the continuum. 
We are now ready to compute 

Z' l/2(^(0)Q(0))o,0 



i+^VK, {pm'u ^^^^^ 

and g/z^ from and flU^ . We have implemented all the series in maple [with 
independent codes in Berlin and Munich] and obtain the corresponding truncated 
series in g^^"^ with coefficients given above and in terms of the moments 

= (0"')o,0, = Ai2. (114) 
Derivatives with respect to Zq are taken with the help of 

d I 



Partial integration implies 

fim+4. + GZQfirn+2-G{m+l)flra = (116) 

which allows us to express all moments by /i = in the final result. In this way 
we have arrived at the coefficients quoted under (jMD, dSSD, dSS]) and (1^ . 

For Zq>0 the integral with ([35]) can be found under 3.323 in [TT] and yields 
the result 

with the modified Bessel function Ki,{.) of index u. For moderately negative Zq 
we have simply summed the expansion for the integrals in z^ to sufficiently high 
order. To compare with ordinary perturbation theory we need the perturbative 
expansion of fi for large Zq, 

11 2 11 34 , , 

/C/Q 

In fact this is a nice pedagogical example for the working of an asymptotic series. 
In Fig. [7] we see the exact /i of f lll7p together with various truncations of the series 
(11181) . There is a range around Zq"^ > 1 where the leading order alone is the most 
decent approximation, a situation reminiscent of the behavior of the perturbative 
series in Fig. |H 
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Figure 7: The behavior of fi from flll7p and its perturbative approximations. 
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